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Abstract 

' Dark strings have recently been suggested to exist in new models of dark matter that 

explain the excessive electronic production in the galaxy. We study the interaction of these 
{Nl ' dark strings with semilocal strings which are solutions of the bosonic sector of the Standard 

bJQ 1 Model in the limit sin 2 9 W = 1, where # w is the Weinberg angle. While embedded Abelian- 
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Higgs strings exist for generic values of the coupling constants, we show that semilocal 
solutions with non-vanishing condensate inside the string core exist only above a critical 
value of the Higgs to gauge boson mass ratio when interacting with dark strings. Above 
this critical value, which is greater than unity, the energy per unit length of the semilocal- 
dark string solutions is always smaller than that of the embedded Abelian-Higgs-dark string 
solutions and we show that Abelian-Higgs-dark strings become unstable above this critical 



value. Different from the non-interacting case, we would thus expect semilocal strings to be 
stable for values of the Higgs to gauge boson mass ratio larger than unity. Moreover, the 
\ one-parameter family of solutions present in the non-interacting case ceases to exist when 

semilocal strings interact with dark strings. 

m ' 
m 

rn 1 1 Introduction 

There is strong observational evidence pQ that approximately 22% of the total energy density 
of the universe is in the form of dark matter. Up until now it is unclear what this dark matter 
should be made of. One of the favourite candidates are Weakly Interacting Massive Particles 
(WIMPs) which arise in extensions of the Standard Model. Recently, new theoretical models 
of the dark matter sector have been proposed [2], in which the Standard Model is coupled 
to the dark sector via an attractive interaction term. These models have been motivated by 
new astrophysical observations [3] which show an excess in electronic production in the galaxy. 
Depending on the experiment, the energy of these excess electrons is between a few GeV and 
a few TeV. One possible explanation for these observations is the annihilation of dark matter 
into electrons. Below the GeV scale, the interaction term in these models is basically of the 
form of a direct coupling between the U(l) field strength tensor of the dark matter sector and 
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the U(l) field strength tensor of electromagnetism. The U(l) symmetry of the dark sector 
has to be spontaneously broken, otherwise a "dark photon" background leading to observable 
consequences would exist. 

Consequently, it has been shown that the dark sector can have string-like solutions, denomi- 
nated "dark strings" and the observational consequences of the interaction of these dark strings 
with the Standard Model have been discussed [I]. 

Topological defects are believed to have formed in the numerous phase transitions in the 
early universe due to the Kibble mechanism [5]. While magnetic monopoles and domain walls, 
which result from the spontaneous symmetry breaking of a spherical and parity symmetry, 
respectively, are catastrophic for the universe since they would overdose it, cosmic strings are 
an acceptable remnant from the early universe. These objects form whenever an axial symmetry 
gets spontaneously broken and (due to topological arguments) are either infinitely long or exist 
in the form of cosmic string loops. Numerical simulations of the evolution of cosmic string 
networks have shown that these reach a scaling solution, i.e. their contribution to the total 
energy density of the universe becomes constant at some stage. The main mechanism that 
allows cosmic string networks to reach this scaling solution is the formation of cosmic string 
loops due to self-intersection and the consequent decay of these loops under the emission of 
gravitational radiation. 

For some time, cosmic strings were believed to be responsible for the structure formation 
in the universe. New Cosmic Microwave background (CMB) data, however, clearly shows that 
the theoretical power spectrum associated to Cosmic strings is in stark contrast to the observed 
power spectrum. However, there has been a recent revival of cosmic strings since it is now 
believed that cosmic strings might be linked to the fundamental strings of string theory [6]. 

While perturbative fundamental strings were excluded to be observable on cosmic scales for 
many reasons [7J, there are now new theories containing extra dimensions, so-called brane world 
model, that allow to lower the fundamental Planck scale down to the TeV scale. This and the 
observation that cosmic strings generically form at the end of inflation in inflationary models 
resulting from String Theory [8] and Super symmetric Grand Unified Theories [9] has boosted 
the interest in comic string solutions again. 

Different field theoretical models describing cosmic strings have been investigated. The U(l) 
Abelian-Higgs model possesses string-like solutions [10 . This is a simple toy model that is 
frequently used to describe cosmic strings. However, the symmetry breaking pattern £7(1) — > 1 
has very likely never occurred in the evolution of the universe. Consequently, more realistic 
models with gauge group SU(2) x U(l) and symmetry breaking SU(2) x U(l) — > U(l) have 
been considered and it has been shown that these models have string-like solutions [11} I12j . 
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Semilocal strings are solutions of a SU(2) g i b a i x U(l)i oca i model which - in fact - corresponds 
to the Standard Model of Particle physics in the limit sin 2 6 W = 1, where 9 W is the Weinberg 
angle. The simplest semilocal string solution is an embedded Abelian-Higgs solution [TT]. A 
detailed analysis of the stability of these embedded solutions has shown [T3] that they are 
unstable (stable) if the Higgs boson mass is larger (smaller) than the gauge boson mass. In 
the case of equality of the two masses, the solutions fulfill a Bogomolny-Prasad-Sommerfield 
(BPS) [2] bound such that their energy per unit length is directly proportional to the winding 
number. Interestingly, it has been observed jT3J that in this BPS limit, a one-parameter family of 
solutions exists: the Goldstone field can form a non-vanishing condensate inside the string core 
and the energy per unit length is independent of the value of this condensate. These solutions 
are also sometimes denominated "skyrmions" and have been related to the zero-mode present 
in the BPS limit. 

In this paper, we consider the interaction of dark strings with string-like solutions of the 
Standard Model in the specific limit sin 2 6yj = 1. The two sectors interact via an attractive 
interaction that couples the two U(l) field strength tensors to each other. This type of interaction 
has been studied before in [15], where the interaction between Abelian-Higgs strings and dark 
strings has been investigated. It has been found that a BPS bound exists that depends on the 
interaction paramater and that Abelian-Higgs strings and dark strings can form bound states. 

Our paper is organized as follows: in Section 2, we give the model, the equations of motion, 
the boundary conditions and the asymptotics. In Section 3, we present our numerical results 
and Section 4 contains our conclusions. 

2 The model 

We study the interaction of a SU (2) 9 ; ^ x U (l)i oca i model, which has semilocal strings solutions 
with the low energy dark sector, which is a U(l) Abelian-Higgs model. 
The matter Lagrangian C m reads: 

(i) 

with the covariant derivatives D^<& = V M <I> — ie\A^<&, = V^ — ie-io*^ and the field strength 
tensors F^ u = d fl A u — d u A^, H^ u = d^a u — d v ap of the two U(l) gauge potential Ap, with 
coupling constants e\ and &i- $ = (0i)^2) T is a complex scalar doublet, while £ is a complex 
scalar field. The gauge fields have masses Mw,i = V^ZiVii * = 1)2, while the Higgs fields have 
masses Mh,% = v / 2Ai^, i = 1,2. The term proportional to e is the interaction term [3]. To be 
compatible with observations, e should be on the order of 1CP 3 . 
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2.1 The Ansatz 



For the matter and gauge fields, we have [TT | [T3 | [TO] : 

<k(p,v)=vihi(p)e iwp , Mp)=mh2( P ) , ap^) = mf(p)e imip (2) 

A u dx>* = — (n - P(p))d<p , a w c^ = — (m - R(p))du> . (3) 
ei e 2 

n and m are integers indexing the vorticity of the two Higgs fields around the z— axis. In the 
following, we will refer to solutions with /i2(p) = as "embedded Abelian-Higgs solutions", 
while solutions with /i2(p) 7^ will be referred to as "semilocal solutions". Note that in the case 
e = 0, the solutions of the semilocal sector of our model are often also denominated "skyrmions" . 



2.2 Equations of motion 

We define the following dimensionless variable x = eiTjxp, which measures the radial distance in 
units of Mp^i/\/2. 

Then, the total Lagrangian C m — ► C m /(r]fe 2 ) depends only on the following dimensionless 
coupling constants 



& = ~2 = 772 2 > 1 = 1 ' 2 ' 9 = — ! 9 = — • (4) 



Varying the action with respect to the matter fields we obtain a system of five non-linear differ- 
ential equations. The Euler-Lagrange equations for the matter field functions read: 

[xH [)' = — + fax{h\ + h\ - l)h x (5) 

x 

{xh <)> = ( n - P ) 2 h2 + ^ + h 2 ^ (6) 

X 

(xfy = — +l32x(f 2 -q 2 )f (7) 

X 

2, fP'\' h?P (n-P)h 2 Rf 2 
(1-e 2 ) — =2^ 2^ J —l + 2eg^— , (8) 

\ X J X X X 

2,fR'\' 2 fR (Ph 2 (n-P)hV\ 
(1-e 2 — =2g 2j — + 2eg[^- y - ^ , 9 

\X J X \ X X J 

where the prime now and in the following denotes the derivative with respect to x. 
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2.3 Energy per unit length and magnetic fields 

The non- vanishing components of the energy-momentum tensor are (we use the notation of [16] ) : 

T£ = e s -e v -e w + u , T Z Z =T$ (10) 

where 

e s = {ti 1 ) 2 + {ti 2 ) 2 + {f) 2 , e v = ^-^ + i-^ 2- , e w = ^^+ 2V ; + ' 



and 



2x 2 2g 2 x 2 g x 2 ' x 2 x 2 ./ 



„ _ ^ ^2 , ^2 ^2 , & ^ 2 ^2 



2 

(11) 



i^ + Z^-lf + f (/ 2 -g 2 ) 2 . (12) 

We define as inertial energy per unit length of a solution describing the interaction of a 
semilocal string with winding n and a dark string with winding m: 

M (n ' m) = J V=g^dxAp (13) 

where is the determinant of the 2 + 1-dimensional space-time given by (t,x,<p). This then 
reads: 

/•oo 

^{n,m) =2n x ( £s + £v+£w+ u ) dx ( U ) 

J 

Note that the string tension T = J \/—gs T^dxdip is equal to the energy per unit length. There 
few special case, in which energy bounds can be given: 

1. For h 2 (x) = 0, the energy per unit length of the solution is given by: 

^(n.m) = 27m r J 2 l9l (p 1 ) + 2nmr, 2 g 2 {f3 2 ) (15) 

where g\ and g 2 are functions that depend only weakly on /?i and (3 2 , respectively. The 
energy bound is fulfilled, when the functions g\ and g 2 become equal to unity. This happens 
at j3\ = (3 2 = 1/(1 — e) and n = m [15] . 

2. For e = and h 2 {x) ^ 0, the energy per unit length of the solution is given by 

/>> m ) = 2nn V l + 2Trmi ] 2 g 2 ((3 2 ) (16) 

where g 2 is a function that depends only weakly on (3 2 with <7 2 (1) = 1. Note that the 
solution of the semilocal sector exists only for f3\ = 1 and fulfills the BPS bound for all 
choices of h 2 (0). 
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The magnetic fields associated to the solutions are given by [15] : 



-P'ix) + £R'(x) , R'(x) 

B z (x) = — 9 —— and b z (x) = , (17) 

e\x e2X 

respectively, where we have used the fact that the part of the Lagrangian containing the field 
strength tensors can be rewritten as [I] : 

- \f^F^ - \h^H^ + ~F„ V H^ => - ~(1 - e 2 )H^H^ (18) 

with G^y = d^A u — d u A^ where = A^ — ea^. The corresponding magnetic fluxes J d 2 x B 
are 

2-7T / e \ / 2-irm 

W = — re m and ip = y 1 — e z , (19) 

ei V 9 J e 2 

respectively. Obviously, these magnetic fluxes are not quantized for generic e. 
2.4 Boundary conditions and asymptotics 

The requirement of regularity at the origin leads to the following boundary conditions: 

hi(0) = , h' 2 (0) = , /(0) = , P(0) = n , R(0) = m (20) 

For /12(C)) = 0, the semilocal strings correspond to embedded Abelian-Higgs strings. Here, we 
are mainly interested in constructing solutions that are truly semilocal, i.e. we require 7*2 (0) 7^ 0. 
The finiteness of the energy per unit length requires: 

/h(oo) = 1 , fta(oo) = , /(oo) = g , P(oo) = , R(po) = . (21) 

The asymptotic behaviour for x — > oo depends crucially on whether the function h,2(x) = 
or /i2(a;) ^ 0. 

1. For /i2(s) = we find: 

P{x oo) = -y^ mi2 [Ci exp (-x/3 + ) + C 2 exp (-x/3_)] + .... (22) 
P(x^oo) = [Ci mn(/3 + )exp(-x/5+) + C 2 mn(/3_)exp(-x/3_)] + ... (23) 

where Ci and C2 are constants, mn(j9±) = (1 — £ 2 )0± — 2 and mi2 = —2eq 2 g. The /?± are 
positive and are given by 



P± - 1 _ £ 2 \ Z V 



The numerical evaluation (see below) shows that for specific values of the coupling con- 
stants the constants C\ and C2 have opposite sign. Hence, the function R(x) can possess 
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a node asymptotically which we have confirmed numerically. However, the numerics has 
shown that these type of solutions exist only for values of e of order one. Hence, we don't 
present them here since we believe that they are unphysical. 

For the scalar fields, we find 

hUx^oo) = 1 + -%exp (-xy/2jh) + — exp(-2x/3+) + — exp(-2x/3_) + ..(25) 

y/X V / X X 

f(x^oo) = q+ -%exp (-xJ2fh) + — exp (-2x/3+) + — exp(-2:r/3_) + ..(.26) 
Jx \ J X X 



C3 and C4 are two constants, while c±, d± depend on the constants C%, . . . , C4 and on (3\ 
and fa- 

2. For h,2(x) 7^ we find : 

2 

P ( x - °°) = ^ + - . ^^°°) = ^2+- (27) 

for the gauge field functions. Here c, cr are constants that depend on the values of the 
coupling constants. For the scalar and Higgs field functions we have 

/ ll (s-oo) = l-^ + ... , h 2 (x) = ^ + ... , /( x -> ro ) = g _^_L_ + .... (28) 

Obviously, the presence of the scalar field /i2(^) changes the asymptotics drastically. While 
for h,2{x) = 0, the gauge and Higgs fields decay exponentially, they have power-law decay for 
h 2 (x)^0. 

2.5 Stability 

Following the investigation in the case e = [13] . we are interested in the stability of the 
embedded Abelian-Higgs string coupled to a dark string. In order to do that we will study 
the normal mode along a very specific (but standard) direction in perturbation space about the 
embedded Abelian-Higgs string coupled to a dark string. We consider the perturbation 

h^x) = h x {x) , h 2 (x) = e^ri(x) , P{x) = P(x) , R{x) = R(x) , f(x) = f(x) (29) 

where the tilded functions denote the profiles of an embedded Abelian-Higgs string coupled to a 
dark string, i.e. solutions to the equations flSJ),©, © and © for fi2(x) = 0. The perturbation 
is denoted by ij and the parameter u> is real. Inserting this perturbation into ([6j) and keeping 
only the linear terms in 77 leads to the linear eigenvalue equation : 

.£ i -±± + V ef ,) V {x)=cJ'T,(z) , V e j f = {n - P J x)? +p 1 (h 1 (xf-l) (30) 



The spectrum of the linear operator entering in (|30p consists of a continuum for u> 2 > and of 
a finite number of bound states (or normalisable solutions) for uj 2 < 0. In the latter case, the 
solutions fulfill 

77(0) = 1 , r/(0) = with tj(x) -» e~^ x for x 00 (31) 

where we have fixed the arbitrary normalisation by choosing 77(0) = 1. 

Only bound states are of interest to us since they signal the presence of an instability. It 
should be pointed out that the functions P(x), h\(x) entering in the effective potential feel the 
effect of the dark sector since the corresponding equations are directly coupled. 

3 Numerical results 

For all our numerical calculations, we have chosen q = g = 1. 

3.1 Stability of the embedded Abelian— Higgs— dark strings 

We have first studied the stability of the embedded Abelian-Higgs strings coupled to dark strings 
by investigating the bound states of (|3U|) for different values of e. Our results for n = m = 1 
and /?2 = 1 are shown in Fig HI 

For £ = 0we recover the result of |13| that the embedded-Abelian-Higgs strings are unstable 
for Pi > 1. For e 7^ 0, we observe that the larger e, the larger the ratio of Higgs to gauge 
boson mass Pi at which the embedded Abelian-Higgs strings coupled to dark strings become 
unstable. In the following, we will denote the value of Pi at which uj 2 = Pf. With view to 
the observations for the e = case, we would thus expect additional solutions with h^ix) 7^ 
for Pi > Pf. In section 3.2, we will discuss the properties of these solutions. 

Let us also remark that our analysis does not reveal the occurence of additional unstable 
modes in the sector explored. 

3.2 Properties of semilocal— dark strings 

In the case e = 0, the two sector do not interact and for the semilocal sector two different types 
of solutions are possible: (a) embedded Abelian-Higgs solutions with /i2(x) = which exist for 
generic choices of Pi |11] and (b) semilocal strings ( "skyrmions" ) with h2(x) 7^ which exist 
only for Pi = 1 [T3]. In the latter case, it was shown that there is a zero mode associated to the 
fact that the energy of the "skyrmions" does not depend on the value of /12(C)). 

The case with e / and h-2,{x) = corresponds hence to the case of an embedded Abelian- 
Higgs string interacting with a dark string. The equations of motion that describe this case are 
exactly those studied in [15]. In |15j . the interaction of a dark string with an Abelian-Higgs 
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Figure 1: We give the value of u 2 (see (f30|) in dependence on (5\ for three different choices of e 
including the non-interacting case e = 0. Here n = m = 1 and 02 = 1- 

string has been studied in detail. Since the only difference between an Abelian-Higgs string and 
an embedded Abelian-Higgs string is the stability - see section 3.1 - we do not discuss this case 
in detail in this paper and focus on the case of semilocal strings interacting with dark strings. 
We have solved the differential equations subject to the boundary conditions numerically using 
the ODE solver COLSYS p2]. 

To see the difference between embedded Abelian-Higgs-dark string solutions and semilocal- 
dark string solutions, we present the energy density Tq , the effective energy density xTq as well 
as the magnetic field B z (see (fT7|) ) in FigJ2]for e = 1/6, (3\ = 3 and /?2 = (l — e)^ 1 = 1.2. Clearly, 
the effective energy density tends to zero very quickly for the embedded-Abelian-Higgs-dark 
string, while for the semilocal-dark string is has a long tail which results from the power-law 
fall off of the functions. Moreover, the magnetic field B z tends to zero exponentially for the 
embedded Abelian-Higgs-dark strings, while it falls off power-like for the semilocal-dark strings. 
Hence, the core of the magnetic flux tube of the latter solution is not well defined. 

While for s = solutions with /i2(x) ^ exist only for (3\ = 1, the situation is different 
here. For e / 0, we find solutions for generic values of /3\, i.e. different from unity. In fact, the 
solutions exist only for (3\ larger than a critical value, /3f r , which depends on the choice of the 
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Figure 2: We give the profiles of the energy density Tq, the effective energy density well 
as the magnetic field B z (see (fTT|) ) for e = 1/6, (3\ = 3 and 02 = (1 — e) _1 = 1-2. We compare 
semilocal-dark string solutions with /i2(0) > (black) and embedded Abelian-Higgs-dark string 
solutions with /12(C)) = (red). 



winding numbers and other coupling constants, in particular e. Moreover, we observe that the 
Pi at which semilocal-dark strings exist is a function of /i2(0). While for e = 0, /?i = 1 for all 
choices of /12(C)), we find that for e / the choice of /12(C)) fixes the value of /3\. 

At the branch of solutions describing a semilocal string in interaction with a dark string 
bifurcates with the branch of solutions describing the interaction of an embedded Abelian- 
Higgs string with a dark string. This is shown in Fig.s [3TH for e = 0.1 and e = 0.5, respectively. 
Note that (3^ is exactly the value at which the embedded Abelian-Higgs-dark strings become 
unstable. 

Here, we give the value of /i2(0) in dependence on (3% for n = m = 1 and P2 = 1-0. Clearly 
at some /3J r , /i2(0) tends to zero which means that h>2{x) = 0. Here the semilocal-dark string 
solutions bifurcate with the embedded Abelian-Higgs-dark string solutions. We also compare 
the energy per unit length of the two types of solutions. Clearly, whenever semilocal-dark string 
solutions exist, they have lower energy than the corresponding embedded Abelian-Higgs-dark 
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Figure 3: The energy per unit length fi^' 1 ^ (in units of 27T7]f) as well as the value of h>2(0) and 
the asymptotic constants c and cr (see (|27p . (|28p ) of the semilocal-dark string solutions are 
shown in dependence on /3i for e = 0.1, 02 = 1 and n = m = 1 (dashed). For comparison, we 
also give the energy per unit length of the embedded Abelian-Higgs-dark string solution (solid). 



string solutions. Moreover, the larger f3i, the bigger is the difference between the two energies 
per unit length. We would thus expect the semilocal solutions to be stable with respect to 
the decay into the embedded Abelian solutions when coupled to dark strings. We also present 
the values of the asymptotic constants c and cr (see (f2"T|h (|2"%|) ). These vanish identically at 

01 = /3f • 

In general, 0f will depend on the choice of 02, n and m: 0f{02, n, m). As shown in [15] in 
the limit h2(x) = a BPS bound exists for /3i = /?2 = (1 — e ) _1 an d n = m. In this limit, the 
energy per unit length (in units of ^nrjf) is just n + m = 2n. We have studied the dependence of 
the energy per unit length on Pi for /3 2 = (1 — e)" 1 where e = 1/6 and e = 0.5, respectively. We 
have chosen n = m = 1. Our results are given in FigEl Interestingly, we find that the branch 
of semilocal-dark string solutions bifurcates with the branch of embedded Abelian-Higgs-dark 
string solutions exactly at 0\= 02 = (1 — E ) ■ For f3\ > (1 — e)" 1 , the energy per unit length of 
the semilocal-dark string solutions is always smaller than that of the corresponding embedded 
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Figure 4: The energy per unit length /i^ 1 ' 1 ' (in units of as well as the value of h 2 (0) and the 

asymptotic constants c and cr (see ([27]) . ([25]) ) of the semilocal-dark string solution are shown 
in dependence on 0i for e = 0.5, 02 = 1 and n = m = 1 (dashed). For comparison, we also give 
the energy per unit length of the embedded Abelian-Higgs-dark string solution (solid). 

Abelian-Higgs-dark string solutions, for 0\ < (1 — e) _1 no semilocal-dark string solutions exist 
at all. Hence, we find that 

0?(p 2 = (l-e)-\l,l) = (l-e)- 1 (32) 

We have also studied the dependence of (3f on the winding of the dark string and the Higgs 
to gauge boson ratio 02 of the U(l) model describing the dark string in more detail. Our results 
are shown in FigJBJ Obviously, 0f increases with increasing s. This is related to the fact that 
the core width of the strings decreases with increasing e. This means more gradient energy 
and hence we have to choose larger values of 0\ to be able to compensate for this increase by 
decrease in potential energy. 

For Pi = 1.0, which in fact corresponds to the BPS limit of the U(l) dark string model for 
e = 0, the value of Pf increases for increasing winding m of the dark string. Again increasing 
m increases gradient energy such that we have to choose larger value of Pi to compensate the 
increase by decrease in potential energy. This is also true when increasing 02- Increasing 02 
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Figure 5: The energy per unit length fjy->^ (in units of ^Trijf) is shown for semilocal strings 
interacting with dark strings as function of f3\ for = (1 — e)^ 1 with e = 0.5 and e = 
1/6, respectively (dashed). For comparison, we also give the energy per unit length of the 
corresponding embedded Abelian-Higgs solutions interacting with dark strings (solid). 

decreases the core size of the dark string, this increases gradient energy and we again have to 
compensate by increasing the value of (3\. 

4 Conclusions 

In this paper we have shown that the interaction of semilocal strings with dark strings has 
important effects on the properties of the former. While embedded Abelian-Higgs strings exist 
for all values of the Higgs to gauge boson ratio when interacting with dark strings, semilocal 
strings with a condensate inside their core exist only above a critical value of the Higgs to gauge 
boson ratio. At this critical value, the embedded Abelian-Higgs-dark strings become unstable. 
The critical value of the ratio depends on the choice of the Higgs to gauge boson ratio of the dark 
string and the windings. In the limit where the ratio tends to the critical ratio, the condensate 
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Figure 6: The value of (5f at which the branch of semilocal solutions bifurcates with the branch 
of embedded Abelian-Higgs solutions is shown as function of e for m = 1, m = 2, respectively 
and /?2 = 1-0, fa = 2.0, respectively. 

vanishes identically and the branch of semilocal-dark string solutions bifurcates with the branch 
of embedded Abelian-Higgs-dark string solutions. Apparently, the presence of the condensate 
lowers the energy in such a way that whenever semilocal-dark strings exist, they are lower in 
energy than their embedded Abelian-Higgs-dark string counterparts. The value of the Higgs to 
gauge boson ratio for which semilocal-dark strings exist depends on the value of the condensate 
on the string axis and increases for increasing values of the condensate. All these results are 
quite different from what is observed in the non-interacting case. In the non-interacting case, 
semilocal strings exist only for Higgs to gauge boson ratio equal to unity and in this limit, the 
energy per unit length is independent of the value of the condensate and in addition fulfills a 
BPS bound. To state it differently : when not interacting with dark strings, semilocal strings 
and embedded Abelian-Higgs strings are degenerate in energy, while the former are lower in 
energy as soon as they interact with dark strings. Since the branch of semilocal-dark string 
solutions bifurcates with the branch of embedded Abelian-Higgs-dark strings at the self-dual 
point of the embedded Abelian-Higgs-dark strings - at which these fulfill an energy bound [15] 
- we expect that semilocal-dark strings are stable. Moreover, they are stable for all choices of 
the Higgs to gauge boson ratio for which they exist and not just - as in the non-interacting 
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case - for Higgs to gauge boson ratio smaller or equal to unity. Since all current observations 
point to the fact that the Higgs boson mass is larger than the gauge boson masses, semilocal 
strings could still be stable when interacting with dark strings. Interestingly - as mentioned 
above - semilocal strings can lower their energy by forming a non-vanishing condensate inside 
their core. This could be important for the evolution of cosmic string networks since next to 
the formation of bound states [18] this would be a further mechanism for the network to loose 
energy. 

We didn't study the gravitational properties of the solutions since we believe that the qual- 
itative features are similar to the case studied in [15]. Since semilocal-dark strings have lower 
energy per unit length than their embedded Ablian-Higgs-dark string counterparts, we would 
expect the deficit angle created by the former to be smaller than that of the latter. Furthermore, 
the critical value of the gravitational coupling at which the solutions become singular is larger 
for the semilocal-dark strings than for the embedded Abelian-Higgs-dark strings. 
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